Triboelectric charging of volcanic ash from the 2011 Grimsvotn eruption 
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Triboelectric charging of different size fractions of a sample of volcanic ash is studied experimen- 
tally. Laboratory experiments demonstrate that the normalised span of the particle size distribution 
plays an important role in the magnitude of charging generated. Previous measurements of the vol- 
canic plumes have shown that ash particles are electrically charged up to hundreds of km away from 
the vent, which indicates the the ash particles continue to be charged in the plume through the 
mechanism of triboelectrification [Harrison et al., Env. Res. Lett. 5 024004 (2010), Hatakeyama J. 
Met. Soc. Japan 27 372 (1949)]. The influence of the normalised span on plume charging suggests 
that all ash plumes are likely to be charged, with implications for remote sensing and plume lifetime. 

PACS numbers: 92.60.Pw, 41.20.Cv 
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Volcanic ash is known to charge electrically, produc- 
ing some of the most spectacular displays of lightning 
on the planet [l|, \^. Although the exact details of ash 
charging processes will vary from one eruption to an- 
other, mechanisms such as triboelectrification, fractoe- 
mission, or the 'dirty thunderstorm' mechanism [l|-|3| are 
all thought to play a role in the electrification of ash near 
the vent. In addition to near- vent charging, observations 
show that charging can also occur in volcanic plumes up 
to hundreds of km from the source region [J, |5| . The sus- 
tained nature of this charge in the presence of electrically 
conducting air, suggests that a self-charging mechanism, 
through the action of triboelectrification, may also play 
a role in the electrification of volcanic ash. Previous the- 
oretical work on triboelectric charging of single-material 
particle systems has shown that the charging is deter- 
mined by the number size distribution [6]. This paper 
details a laboratory investigation into triboelectric charg- 
ing of a sample of ash from the Grimsvotn eruption in 
Iceland in 2011, in terms of the particle size distribution, 
using specially designed apparatus. 

Charging arising from contact between two different 
material surfaces can be understood as a result of the 
different work functions of the materials, however tribo- 
electric charging in systems of identical materials cannot 
be explained in this way. Lowell and Truscott presented 
a model for triboelectric charging between macroscopic 
samples of identical materials based on spatial localisa- 
tion of electrons on the material surface [7] . Spatial lo- 
calisation of electrons prevents relaxation of electrons in 
high energy states to vacant low energy states elsewhere 
in the material. Contact between two surfaces provides 
a relaxation mechanism where a localised high energy 
electron on one surface can move to a vacant low energy 
state on the other surface, resulting in electron transfer 
between surfaces. 

This model has more recently been developed to de- 



scribe triboelectric charging of granular systems |6] . The 
number of trapped high energy electrons is assumed to be 
proportional to the particle's surface area, i.e. the sur- 
face charge density is the same for all particles, and the 
number of low energy electrons is zero. In a collision, a 
high energy electron in one particle will be transferred to 
a low energy state in the other particle. If both particles 
have equal numbers of high energy electrons, there is no 
net charge transfer. However, if only one particle has a 
high energy electron, this will be lost to the other parti- 
cle. Smaller particles will therefore lose all their trapped 
high energy electrons before the larger particles, while 
continuing to receive electrons into vacant low energy 
states, causing net electron transfer from large to small 
particles. This results in an average negative charge on 
the smaller particles and average positive charge on the 
larger particles. 

Lacks and Levandovsky present simulations of parti- 
cle dynamics to illustrate their model, which reproduces 
the negative charging of smaller particles and positive 
charging of larger particles observed empirically [y| . This 
model has since been extended to include geometric con- 
siderations which favour electron tunnelling from large 
particles to small [8|. In addition to these numerical 
studies. Forward et al. present experimental studies of 
triboelectric charging in soda lime glass. Mars and lu- 
nar regolith simulants and demonstrate that the particle 
size dependent charging seen in natural phenomena (e.g. 
dust devils, volcanic plumes) can be reproduced in the 
laboratory |9l4ll|. 

Pahtz et al. present a model of electron transfer be- 
tween identical dielectric grains in an electric field |12| . 
The applied field polarises the grains and when two oppo- 
sitely charged surfaces collide, electrons are transferred. 
Following separation, the applied field repolarises the 
grains. This charging model is not applicable to our 
experiments as care is taken to ensure there are no ap- 
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FIG. 2: Volumetric particle size distributions measured with 
the Malvern Mastersizer 2000. (a)-(c) sieved samples, with 
the size fractions defined by the sieves indicated (d) artifi- 
cial narrow bimodal distribution (50:50 mixture of 45-63 /im 
and 90-f 25 /ini) (e) artificial wider bimodal distribution (50:50 
mixture of 45-63 /im and f25-f80/im) and (f) the size distri- 
bution of the sample before sieving. 



FIG. 1: Scanning electron microscope image of ash sample. 
The non-spherical nature of the particles is clear. Image cour- 
tesy of David Pyle. 



plied external electric fields, however it may contribute 
to charging in plumes. 

The 2011 eruption of the Grimsvotn volcano began on 
21 May, and the eruption was associated with consider- 
able volcanic lightning [13]. Preliminary estimates have 
put the total amount of tephra ejected at 0.25 km^ |14| . 
Ash was collected on 26 May at Kirkjub^jarklaustur, 
approximately 75 km SSW of the Grimsvotn crater. The 
eruption ended on 28 May. Scanning electron microscope 
images of the sample (Figure [T]) show the particles to be 
angular and to have a wide range of sizes. 

Ash diameter distributions were measured with a 
Malvern Mastersizer 2000, which uses laser diffraction to 
calculate volumetric size distributions of suspended sam- 
ples. Volumetric size distributions were obtained for the 
sample before and after dry sieving (to preserve aggre- 
gates that could contribute to the plume's electrostatic 
properties) to separate different size fractions, shown in 
Figure [5] and summarised in Table HI Dry sieving showed 
that the larger ash particles were slightly darker in colour 
than the smaller particles, which suggests the ash sample 
is made of a mixture of different substances. These dif- 
ferent substances may triboelectrically interact with each 
other, in addition to the charge transfer as a function of 
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Before sieving 


1.919 


0.832 


45-63 ^m 


1.801 


0.859 


63-90 ^m 


1.642 


0.872 


90-125 /im 


0.852 


0.906 


125-180 fim 


0.775 


0.909 


45-63 and 90-125 ^m 


1.144 


0.892 


45-63 and 125-180 fim 


1.469 


0.877 



TABLE I: Volumetric size distribution summaries. The nor- 
malised span is a non-dimensional index of the polydisper- 
sity of the distribution, defined by the normalised interdecile 
range |15l |. The modality coefficient b is calculated from the 
skewness and kurtosis of the distribution. A value of b greater 
than 0.55 indicates a multimodal distribution ^16] . 



size. 

Three sieved samples were tested, from nominally 45- 
63, 63-90 and 90-125 /im distributions, defined by the 
sieves. Two artificial size distributions were created from 
mixing 50:50 samples (by mass) of 45-63 /iin and 90- 
125 /im, and 45-63 /im and 125-180 /iin, to generate a nar- 
row bimodal and a wide bimodal distribution. 

Here we use the term normalised span, a non- 
dimensional index of the polydispersity of the distri- 
bution, defined as the difference between the 90*'' and 
10*'' diameter percentiles, divided by the median diam- 
eter [15|. As might be expected, the pre-sieve size dis- 
tribution showed the greatest polydispersity, with most 
of the particles between 20 and 200 /im. The sieve and 
Malvern sizer diameters only agree approximately, which 
may be due to the assumption of sphericity used by the 



Malvern instrument; this is clearly incorrect as demon- 
strated by Figure [Ij The sieved samples size distribu- 
tions all show a substantial tail of fine particles, which 
we believe cannot be from cross-contamination within the 
particle sizer, due to careful experimental technique. An 
alternative explanation for the fine tail could be disag- 
gregation whilst the samples were in the Malvern sizer, 
however, a similar fine tail would be expected across all 
the measurements, which was not seen (Figure [2]). To 
quantify the modality of the samples, the modality coef- 
ficient b was calculated using the skewness and kurtosis 
of the distribution. A value of b greater than 0.55 indi- 
cates a bimodal or multimodal distribution [16] . Using 
the modality coefficient and span allows the samples to 
be divided into three types: bimodal with a narrow span, 
bimodal with a broad span and monomodal with a broad 
span. 

Electrostatic charging of the Gri'msvotn ash was in- 
vestigated using a grounded tube, through which ash is 
dropped onto an isolated Faraday cup connected to a sen- 
sitive electrometer, shown in Figure [3l In an optimised 
delivery technique, ash was delivered to the charge appa- 
ratus via aluminium tubes with a sliding floor mounted 
on a rotating metal turntable, supported vertically above 
the inlet. The charge, AQ, transferred from the ash to 
the Faraday cup is related to the change in voltage, AV 
measured at the cup and the capacitance, C, of the sys- 
tem (130 pF), by the relationship AQ = CAV. The 
change in Faraday cup voltage was recorded by a Camp- 
beh CR3000 data logger at 300 Hz. 

For each ash charging experiment, 0.2 g of Gri'msvotn 
ash (baked to removed water) was weighed and trans- 
ferred to an individual delivery tube. To minimise the 
unwanted effect of self-charging of the ash during han- 
dling, ash was left in the delivery tube for 30 minutes 
before each test, to allow time for any residual charge 
to decay. We therefore assume that any charge mea- 
sured on the ash after descent is entirely from triboelec- 
tric charging during interactions whilst the ash is falling 
under gravity, analogously to a volcanic plume in the at- 
mosphere. We did not observe the any aggregates in the 
samples before or after the ash drop experiments. 

Charging experiments were undertaken with the flve 
different size distributions of Grimsvotn ash described 
above: the three sieved fractions and the two artificial 
distributions. For each size distribution, ten ash charg- 
ing experiments were performed to improve the sampling 
error. 

Figure|4]shows a typical charging trace measured at the 
Faraday cup, with the sign of the charging inverted for 
clarity. The ash is released a.t t — 0. Initially the charge 
decreases to a minimum value before increasing exponen- 
tially to a maximum value. The change in charge mea- 
sured by the Faraday cup for the three types of samples 
are summarized as box and whisker plots in Figure [SI and 
their distributions compared using the Wilcoxon test. 
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FIG. 3: Schematic diagram of the ash charge apparatus show- 
ing the ash delivery apparatus at the top, and collecting Fara- 
day cup at the base. 



Figure [5] shows the net charge difference for all three 
groups. For the two bimodal distributions, the nar- 
row span group shows much smaller charging (median 
value 13.3 pC) than the broad span group (median value 
28.7 pC) to better than 95% confidence. Including the 
broad monomodal distribution in the broad span group 
does not alter this result. 

For the samples with broad span, there is no difference, 
again to better than 95 % confidence, between samples 
with a monomodal or bimodal distribution where the 
median values of the charge difference are 41.6 pC and 
28.7 pC. This demonstrates that the normalised span of 
the distribution may affect the magnitude of the charg- 
ing, while the modality of the sample does not. The 
results obtained by comparing samples grouped in terms 
of span and modality rather than size also indicates that 
the particle distribution has a greater effect on charging 
than any effects of varying composition with size. 

In conclusion, charging experiments show that 
Grimsvotn ash is easily electrified via the self-charging 
mechanism, with the span of the particle size distribution 
playing an important role in the magnitude of the charge 
generated. Samples with the largest normalised span in 
particle sizes {i.e. a variety of different sized particles) 
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FIG. 4: A typical Faraday cup charge measurement trace with 
the sign of the charging inverted for clarity. Three points are 
marked and the charge differences between these points are 
used to compare the different experiments. The net charge 
differences are summarised in the box plot shown in Figure [5] 
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FIG. 5: Net charge difference at the Faraday cup for the three 
samples. The central mark in each box (red, color online) 
shows the median change for each group. The edges of the box 
show the 25"" and TS*** percentiles and the whiskers extend to 
1.5 times the inter-quartile range. The data points are shown 
as black points and data points outside this limit are shown 
as (red) crosses. 



were observed to generate the largest magnitude charges. 
It is also observed that the span of the particle distribu- 
tion dominates over the modality (either monomodal or 
bimodal) . These results support the theory of triboelec- 
tric charging in granular systems proposed by Lacks and 
Levandovsky [G]. 

These findings have implications for the remote sensing 
of volcanic ash via electrostatic techniques as the amount 
of charging will change with the particle size distribution, 
giving different charging behaviour in different eruptions, 
in different phases of an eruption and as the particle size 
distribution changes through gravitational settling. Sus- 
tained triboelectric self-charging of volcanic plumes dis- 



tant from the vent is possible as long as there is a distinct 
particle size distribution. Further work is needed to de- 
termine the minimum span needed for sustained charg- 
ing. 
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